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Abstract 

In the Standard Model, the box-diagram mediated decay b — * ssd is predicted 
to occur with an extremely small branching ratio of below 10 -11 , thus providing a 
safe testing ground for exposing new physics. We study this process in several Two 
Higgs Doublet Models and explore their parameter space, indicating where this process 
becomes observable. 
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In the study of possible virtual effects in B decays which are due to physics beyond the 
Standard Model(SM), a crucial task is to subtract reliably the SM contributions. For the 
best studied rare b decay process, b — > s + 7, many efforts have been made in the last 
decade to calculate it as accurately as possible within SM. However, since the signals of new 
physics in b decays emerge unluckily with branching ratios at the level of 10~ 7 or less the 
commonly focused processes of which b — > s + 7 is the prototype are very problematic for 
testing new physics against the SM background. Moreover, other yet unobserved B decays, 
like various B — > r processes, have also been shown to be rather insensitive to a large class 
of new physics models @. 

In a recent letter [[J, we have emphasized an alternative approach to the challenge 
of identifying virtual effects from new physics in b decays, by the consideration of processes 
which have negligible strength in the SM. Such processes could thus serve as sensitive probes 
for new physics, relatively free of SM "pollution" . In we focused on the b — > ssd transition, 
which is a box-diagram induced process in the SM with a very small branching ratio of below 
I0~ n , and we studied possible effects from the minimal supersymmetric standard model 
(MSSM) and from a supersymmetric model with R-parity violation. We have shown that 
the existing limits on parameters of MSSM and of R-parity violating interactions allow this 
process to occur well in excess of the SM rate, thus providing an unusual opportunity for 
stricter limits or, hopefully, for discovering new physics effects. In the present letter we 
undertake a study of the occurrence of the b — > ssd transition in Two Higgs Doublet Models 
(THDMs)@, ^, ||, which are frequently considered as likely candidates for the extension of 
the SM. We shall study two models in which the charged Higgs exchange is contributing to 
box diagrams, which we denote as Model I [|J and Model II f5|, respectively, as well as a 
third model allowing a tree level transition mediated by neutral Higgs bosons 0. 

We begin with some comments on the calculation of the W-box diagrams in the SM. 
Due to strong cancellations between the contributions from the top, the charm and the up 



2 



quarks in the loops, the leading order SM result for the b — > ssd decay rate is 
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with x = m^/m 2 , y = ml/m^. The first term in ([J) is suppressed by the small Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements VtdV t * s , while the second term is suppressed by 
y = m\lm\y and contributes about a half of the CKM suppressed term. In this second con- 
tribution, we have neglected a kinematics dependent term when we perform the integral over 
the loop momentum. This amounts to neglecting a small (ml/m 2 v ) \n(f(p)/m c ) contribution 
with f(p) a function depending on the external momenta p. We have checked the effect of 
the neglected dependence numerically and found that it never exceeds 10% of the ml/m^r 
term in the whole kinematic regionQ. Since the involved CKM matrix elements are not well 
bounded and the relative phase of the two terms is not fixed, we can only determine a range 
for the branching ratio of b — ► ssd, which turns out to be always below 10~ n in the SM. 
Note that QCD corrections may not change the value by orders of magnitudes, if compared 
with the analogous processes B°B° and K°K° mixing 0. All these features combine to 
single out this process as a very sensitive one to new physics. 

The THDMs are the minimal extensions of the SM. With one more Higgs doublet, 
one has to suppress the tree level flavor changing neutral current (FCNC) interactions due 
to neutral Higgs bosons to be consistent with the data. The Model I allows only one Higgs 
doublet to couple to both the up- and the down-type quarks ||. In the Model II one Higgs 
doublet is coupled only to up-type quarks while the other doublet is coupled only to down- 
type quarks |J, and thus the Higgs content in the Model II is the same as in the MSSM. In 
2 We note that this dependence does not appear in the calculation of the K K mixing, where the external 
momentum squares are of the order of m 2 and can be safely neglected compared to m 2 , nor in the case of 
BB mixing, where even the contributions of the order of m 2 /??! 2 ^ are sub-dominant. 
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Figure 1: The ratio R of the amplitude induced by the rnl/m^r term over that by the CKM 
suppressed one. The solid line, dashed line, dash-dotted line and dotted line correspond to 
m H + = 100, 200, 400, 800 GeV, respectively. 

both the Model I and Model II, discrete symmetries are enforced to forbid the more general 
couplings and thus the tree level FCNC interactions are absent. In another model, the Model 
III H, the tree level FCNC is assumed to be small enough to lie within the experimental 
bounds. 
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The relevant Lagrangian in the Model I and II is 

£ = ^^[cot(3u i>R (M u V) ij d j}L -Cu i , L (M D V) ij d jjR }H + + h.c, (4) 

V 2 My/ 

where V represents the 3x3 unitary CKM matrix, My and Mq denote the diagonalized 
quark mass matrices, the subscripts L and R denote left-handed and right-handed quarks, 
and i,j = 1,2,3 are the generation indices. For model I, £ = cot/3; while for Model II, 
£ = — tan (3. 
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In the charged Higgs mediated box diagrams, the couplings proportional to the up 
quark masses dominate the contribution to b — > ssd decay. In both models we have 
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with z = m 2 H+ /m 2 . The functions A(x,z) and B(x,z) denote the new contributions from 
the charged Higgs. Because we have dropped in (^j) the terms which are proportional to a 
factor less than morris /M^ at the amplitude level, the limit cot (3 — > corresponds to the SM 
case. As in the SM, the result depends on the relative CKM phase between the terms with 
the m 2 /m 2 v and the CKM suppressed contributions. In Fig. 1 we plot R as the function of 
cot (3, where 



R 



\yV cd V£[g(x,y) + B(x,z) 



\V td V t * s [f(x)+A(x,z)]\ 
is the ratio between the amplitudes induced by the m^/m 2 ^ and by the CKM suppressed 
terms. It is clear that the SM is the limit with the maximal importance of the m^/m 2 ^ 
contribution, while in the limit of large cot (3 this contribution is negligible. 

In the numerical calculation, we use m& = 4.5 GeV, m c = 1.5 GeV, \V ts \ = 0.04, 
\V t d\ = 0.01 and \V cd \ = 0.22. We set the relative phase between the two terms in (|) as 
zero which corresponds to a maximal constructive interference. In Fig. 2, we show the 
numerical results for branching ratio of b — > ssd decay as a function of cot (3 for different 
charged Higgs masses. We have taken m#+ as 100-800 GeV while cot (3 ranging from 0.1 to 
5. Some semi-quantitative considerations suggest a small cot (3 (< 2)|] in which region the 



We refer the reader to H for more detailed discussions. 
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Figure 2: The branching ratio of b — > ssd as a function of cot (3. The solid line, dashed line, 
dash-dotted line and dotted line correspond to m H + = 100,200,400,800 GeV, respectively. 
Regions above the lines are excluded. 

decay branching ratio for b — > ssd decay is unobservable. However, phenomenologically a 
cot (3 as large as 5 is still consistent with the low energy data and the direct search of the 
Higgs boson in top decays. In this parameter space of a large cot (3 the charged Higgs box 
diagrams are sizable. We note that in the MSSM cot (3 < 1 is preferred in model building, 
corresponding to an unobservable contribution from the charged Higgs box diagrams in this 
model §. 

As can be seen in Fig. the branching ratio in the Model I and II can be of the order 
10~ 8 at large cot (3 region. The searching of the decay b — > ssd and its hadronic channels 
— > K ± K ± X in B experiments will further constrain the parameters in the THDMs. 
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In the THDM III, there exists the Yukawa interaction |J 

C = £,ijQi,L<t>2Dj t R + (the up quark sector) + h.c. (9) 
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Here, the scalar Higgs doublet 02 mediates the FCNC transitions di «-> dj at the tree level, 
if the coupling £y is nonzero. As discussed in the literature, FCNC effects induced by the 
loop diagrams are always negligible compared to the tree level ones in the Model III [pL0 | , 



hence the box diagrams with charged Higgs can be safely dropped in the present decay. We 
will neglect the unimportant QCD correction since no GIM-like cancellation happens in the 
process. 

The couplings in (Q) are constrained strongly from the neutral meson mixing by the 
requirement that the FCNC contribution to the mixing from the interactions in (|9|) does not 
exceed its experimental value. In the KK and B S B S mixing, the dominant contributions are 
proportional to ^ d and respectively |10 



* - K (f + f ) , do) 



with F = K, B s , and q = d,b, respectively, and 



M§ = — I fpMp + 



6 y (m s + m q ) 2 J 

Here mu and vtia are the masses of the neutral scalar and pseudoscalar Higgs bosons. Note 
that we have taken £y = ^ while it is not difficult to generalize to the case without this 
assumption. In numerical calculations, we use fx = 160 MeV, fs a = 200 MeV, AMk = 
3.491 x 10~ 15 GeV [11]. Taking itla = rah = rnn, we obtain the bound 



— < 8.3 x lO^GeV" 1 (12) 
m H 

from AMk- The experimental lower limit from B s — B s mixing 

AM Bs > 5.2 x 10~ 12 GeV [12], (13) 

gives no constraint on £, s b/mH, since the contribution in the SM itself exceeds this number. 
Furthermore, free from any assumption made for the FCNC Higgs couplings in the lepton 
sector, the bounds from B s — ► klj and B —> Kklj (k,lj = e,/i or r) do not exclude even 
Csb/nT-H as large as 10 _1 fL3 |. 
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Figure 3: The branching ratio of b — > ssd in THDM III as a function of \^bs^sd\/ m< H- 

In the presence of the interactions @, b — > ssd can be induced by a tree diagram 
exchanging the neutral Higgs bosons h (scalar) and A (pseudo-scalar), with the amplitude 



A = % - isbisd \\{sb){sd) \{s^b){s^ h dy\ . 
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The decay rate is thus 

r = 
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The numerical results are shown in FigfJ as a function of \Csb^sd\/fn 2 H- 

From Fig. 3, it can be seen that the decay b — > ssd is observable in the Model III, 
if \£,sb£,sd\ / rn 2 H > 1CT 10 GeV -2 , which corresponds to a branching ratio around 10~ 9 . This 
requires at least > 10~ 3 GeV -1 . Corresponding to this number, the neutral Higgs 

contribution to AMb s is 10 6 times larger than its present lower limit. To our knowledge, 
such a large AM Bs is difficult to exclude. 
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In summary, we have presented a study of the b — > ssd process in several THDMs. Firstly, 
we confirmed that the charged Higgs box contribution in MSSM is indeed negligible H, 
while on the other hand in Models I 0] and II || this contribution can induce observable 
effects at the 10~ 8 level for the branching ratio. A large B S B S mixing, which is at least 
10 6 times larger than its present lower limit, is required in Model 3 for this process to 
be observable. Combining the above results with those of Ref ||, we reemphasize that the 
search for the processes we recommended ||, like B~ — > K~K~n + , will serve immediately 
to get better limits on the parameters of R-parity violating theories; at the later stage, when 
lower experimental limits are obtainable, it would constitute a direct check of the various 
Non-Standard Model theories investigated here and in Ref [BJ. 

The work of KH and DXZ is partially supported by the Academy of Finland (no. 
37599). The work of PS is partially supported by the Fund for Promotion of Research at 
the Technion. 
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